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ABSTRACT 

The Unified Model of Active Galactic Nuclei (AGN) has provided a successful expla- 
nation for the observed diversity of AGN in the local Universe. However, recent analysis 
of multi-wavelength spectral and image data suggests that the Unified Model is only a 
partial theory of AGN, and may need to be augmented to remain consistent with all 
observations. Recent studies using high spatial resolution ground- and space-based ob- 
servations of local AGN show that Seyfert class and the "core" (r < 1 kpc) host-galaxy 
morphology are correlated. Currently, this relationship has only been established quali- 
tatively, by vi sual inspection of the core rnorph ologies of low redshift {z < 0.035) Seyfert 



host galaxies (jMalkan. Gorjian and Tamlll998l ). We re-establish this empirical relation 



ship in Hubble Space Telescope (HST) optical imaging by visual inspection of a catalog 
of 85 local {D < 63Mpc) Seyfert galaxies. We also attempt to re-establish the core 
morphology-Seyfert class relationship using an automated, non-parametric technique 
that combines both existing classification parameters methods (the adapted CAS, G- 
M20), and a new method which implements the Source Extractor software for feature 
detection in unsharp-mask images. This new method is designed explicitly to detect 
dust features in the images. We use our automated approach to classify the morphology 
of the AGN cores and determine that Sy2 galaxies visually appear, on average, to have 
more dust features than Syl. With the exception of this "dustiness" however, we do 
not measure a strong correlation between the dust morphology and the Seyfert class 
of the host galaxy using quantitative techniques. We discuss the implications of these 
results in the context of the Unified Model. 



1. Introduction 

Active Galactic Nuclei (AGN) are sustained by the accretion of material from their local 
environment onto a super-massive (M > 1O^-1O''M0) black hole. In the Unified Model of AGN, 
the observed diversity in emission-line profiles of AGN is believed to be an observational bias 
introduced by the relative inclinations (with respect to the observer) of the central engiiie as it 



is nested within a toroid of dense molecular material (jBarthel et al.l Il984l : lAntonucci et al.l Il993l ) . 
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Observations of the "zoo" of AGN (e.g., Seyferts, BL LAC objects, Radio galaxies) from X-ray 
to r a.dio wavelengths have b een remarkably well-explained by the Unified Model (for a review, see 



e.g., 



Urrv &: Padovanilll995l ). 



Despite the success of the model, numerous AGN in the local Universe are not well-explained 
within the paradigm of the Unified Model. Many tests of the Unified Model have concentrated on the 
observed diversity in the properties of Seyfert galaxies, which are broadly classified by their emission 
line profiles as: a) Syl-1.9 (Syl), observed with both broad {v > lO^km s~^) and narrow line 



emission; and b) Seyfert 2 (Sy2), observed only with narrow line emission. For example. iTran I (|2001 



20031 ) identified Sy2 AGN that lack "hidden" Syl AGN as predicted b y Unified Model, ind i catin 



that Sy2s may not be — as a class — identical to Syl AGN. Furthermore. iPanessa fc Bassanil (j2002l ) 



found that the column density of absorbers in Sy2 AGN implies the exist ence of dust absorbers on 



a larger physical scale (r > 1 kpc) than the molecular toroid. Recently, iRicci et al.l (120111 ) found 
"excess" X-ray emission from reflection in Sy2 AGN, that did not appear to a comparable extent 
in Syl AGN, indicating an environmental distinction between these two classes of AGN. 

Malkan, Gorjian, & Tam (1998, hereafter MGT98) tested the Unified Model via a "snapshot" 
campaign (see ij2] for details) conducted with HST Wide Field Planetary Camera 2 (WFPC2) 
using the F606W (Ao=5907A) filter, in which they observed the morphology of the inner core (~1 
kpc) of 184 local (z < 0.035) Seyfert galaxies. The authors visually inspected these images and 
determined that Syls are preferentially located in galaxies of "earlier-type" core morphology, and 
conversely that Sy2 AGN are more often hosted by galaxies with "later-type" cores. MGT98 also 
determined that the distribution of dust is more irregular and extends closer to the nucleus in 
Sy2 galaxies than it does in Syl AGN. Hereafter, we refer to these two empirical relationships as 
the "MGT98 relationship." These independent studies suggest that there may be a fundamental 
physical distinction between Syl and Sy2 galaxies, that is not explained by the relative inclination of 
the thick, gas-rich toroidal in which the AGN central engine is embedded. The contemporary debate 
on the nature of AGN is not framed exclusively by the Unified Model; other model s of the central 



engin e and the dusty accretion disk do exist (e.g., the "clumpy torus" model of iNenkova et al 



20081 ). but we will discuss our analysis in the context of the Unified Model to provide an easier 



comparison with published results in the literature. 

In the past two decades, HST images and spectroscopy has been used to study hundreds of low 
redshift (z < 0.1) AGN, and many of these observations are now available in the Hubble Legacy 
Archive (hlaH an online repository maintained as a service for the community by the Space 
Telescope Science Institute. 

In this study, we test the Unified Model using images downloaded from the HLA. Specifically, 
we re-examine and extend the analysis first established in MGT98 using a catalog of 85 Seyfert 
galaxies selected using the criteria outlined in ^ In S}3l we present the results of the visual in- 
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spection and classification of our catalog Seyfert galaxies. In ^ we present, apply and discuss an 
automated technique, which we use to quantify the distribution of any dust features (e.g., dust, 
stellar clusters, etc.) present in the cores of our catalog galaxies. This classification technique 
quantifies the distribution of the dust features that were used to qualify the degree of dust irregu- 
larity or morphological class of a galaxy in the original visual inspection in ^ In ^ we present 
a new automated technique developed to detect the dust features, which were identified in ^and 
used in the visual classification of the galaxies' cores. We discuss the results, and implications, of 
our qualitative visual and quantitative automated analysis in ^ T hroughout, we assume a ACDM 
cosmology with ^7^=0.27, r2A=0.73, and Ho=70 km"^ s"^ Mpc"^ (jKomatsu et ahlboill ). 



2. Data and Image Processing 

To test the MGT98 relationship, we require a sufficiently large sample of Syl and Sy2 AGN 
to ensure that any result can be interpreted in a statistically meaningful way. We therefore use the 
following selection criteria to identify this sample of AGN: 



• Initial Catalog: We develop a large (N~240) catalog from three large HST surveys of Seyfert 
galaxies (Ho et al. 1997; MGT98; Ho & Peng 2001) that were included in the NASA/IPAC 
Extragalactic Database (neeEI). We refer the reader to the respective surveys for specific 
details associated with the sample selection of these AGN. Together, these surveys can be 
used to produce a catalog that is generally representative of the morphological diversity of 
Seyfert galaxies, although none of the samples is strictly volume complete. 

• HST WFPC2 F606W HLA Images: At optical wavelengths, the resolution of features with 
a linear spatial extent of 10 < r (pc) < 100 can only be achieved with large-aperture space- 
based observatories. Thus, we required our galaxies to have HST WFPC2 F606W filter 
images in the HLA, prepared as mosaics of the WFl-3 and PC CCD images multidrizzlecji] to 
a uniform O'.'IO pixel scale. We only used the mosaiced images to ensure that the intrinsically 
different pixel scales of the individual CCDs did not bias the identification and classification 
of sources. The HLA contains an image for more than 90% of the galaxies included in our 
initial catalog with this specific camera and filter combination. The fact that these images 
are available is partly a selection bias. Many of the observations we include in our catalog 
were observed by MGT98 in the snapshot campaign. Note that the F606W filter samples 
longward of the 4000A break at all relevant redshifts in our catalog. This broad filter includes 
the rest-frame Ha and [Nil] line emission which, in AGN, can be prominent. In ^we discuss 
the effect of this emission on our qualitative analysis. 
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• "Face-On": We only included "face-on" galaxies to ensure that the dust features classified 
in ^ are physically confined to a region relatively close to the core (1 kpc) of the galaxy. 
We estimated the angle of inclination by eye, and excluded an additional 20% of AGN that 
appeared at inclinations approximately greater than 30°. We did not exclude those galaxies 
with inclination angles that could not be estimated (i.e., irregular galaxies), nor do we exclude 
elliptical galaxies. 

• Distance less than 63 Mpc: We are interested in characterizing the structural properties of dust 
features with a linear size scale greater than 100 pc (for more details, see ^5.ip . We require at 
least 3.5 WFPC2 pixels (0'.'35 in the HLA mosaic images) to span this physical scale. This sets 
the maximum allowable distance to a catalog galaxy of 63 Mpc or, equivalently, to a redshift 
z < 0.015. Sub-kiloparsec scale features (e.g., dust lanes and clump-cloud formations such 
as bars, wisps, and tidal features like warps and tails) are easily discernible in galaxies nearer 
than this distance observed at the HST spatial resolution. We excluded an additional 50% of 
galaxies that were at distances greater than 63 Mpc. We model and discuss the dependence 
of the morphological classification parameters on spatial resolution (Appendix A), and the 
galaxy distance (Appendix B). 



31 Syl and 54 Sy2 galaxies from the initial sample met all of these selection criteria, combined 
for a total of 85 Seyfert galaxies. This large sample ensures that the (Poisson) uncertainties from 
small number statistics are small. Our catalog includes significantly fewer Syl than Sy2 galaxies, 
partly due to a bias towards Sy2 AGN in the initial sample. For example, only 44% of the galaxies 
in MGT98 are classified as Syl AGN. In the Unified Model, this represents a bias in the opening 
angle through which the AGN is viewed. Though this bias may be present, it will not significantly 
affect this study, because we are investigating the core morphological distinctions between the AGN 
sub-classes of the host galaxies (i.e., on scales of hundreds of parsecs, well beyond the ~parsec scale 
of the thick, dusty torus). Where the data are available from NED, we plot the number of Seyfert 
galaxies by their deVaucouleurs galaxy type and 60/im fiux in Figures [T] & [2l respectively. These 
figures demonstrate that our catalog is not strongly dominated by a particular galaxy type or 
observed AGN luminosity. 

We prepared the HLA mosaiced images for analysis by first visually identifying the (bright) 
center of each galaxy. We extracted a core region with physical dimensions of 2x2 kpc centered 
at this point. The HLA images that we used have only been processed to the Level 2 standard, 
i.e., only images acquired during the same visit are drizzled and mosaiced in the HLA. Many of 
the galaxies were originally imaged as part of HST snapshot surveys (single exposures with texp — 
500s). As a result, cosmi c rays can be a sign ificant source of image noise in the mosaics. We used 
the routine, l.a.cosniic ( van Dokkum 2001 ) to clean the CCD images of cosmic ray^. Initially, 



we implemented l.a.cosmic using the author's suggested parameters, but found by iteration that 
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a lower value for the object— detection contrast parameter, sigclip=2.5, produced cleaner images 
without significantly affecting the pixels of apparent scientific interest. Additional cleaning and 
preparation of the imaging was necessary for the following analyses, and we discuss those task- 
specific steps taken in ^4.11 

3. Visual Classification of Core Morphology 

The core morphologies of the AGN-host galaxies are diverse and early- and late-type mor- 
phologies, with varying degrees of complexity in dust and gas features, are represented in the 
catalog. In Figure[3l we provide images of a subset (12) of our galaxies; each image has been scaled 
logarithmically. Images of all (85) galaxies are available in the electronic supplement to this article. 
Here, we use this subset of galaxies specifically to discuss the various dust features and structures 
that we classify by eye. 

Galaxies in our catalog display a wide variety of spiral arms-like features. In Figure El we 
provide images of two galaxies (MARK1330, NGC3081) that show distinct spiral arms. These 
features are easily distinguished from the ambient stellar light profile in the core of each image. In 
some galaxies, these arms are reminiscent of galactic-scale spiral features, such as barred spirals 
(NGC3081). Some spiral-arm like features are more unique. For example, MARK1330 has a single 
arm that appears to in-spiral and connect with the bright core of the galaxy. Furthermore, some 
galaxies appear to be relatively dusty with numerous features of various size scales, appearing in 
either organized or chaotic features (e.g. NGC1068, NGC1386, NGC1672, NGC3393). 

In Figure EJ we also provide examples of galaxies whose cores are relatively sparsely populated 
with dust features. In some cases (e.g., NGC3608), these galaxies have few dust features. In other 
galaxies (e.g., NGC1058) dust features appear most pronounced in the core of the galaxy (100-200 
pc) and are less significant at large radii. 

We have visually inspected and classified each of the 85 galaxies in our catalog, first using the 
following criteria that were defined and used in MGT98. We divide these criteria into two general 
classes: 

Class 1 — Dust Classifiers: 

• DI : Irregular dust; 

• DC : Dust-disk/Dust-lane passing close or through center (i.e., bi-sected nucleus); 

• D : Direction of dust lanes on one side of major axis, where direction is N, S, E, W, NW, 
NE, SW, or SE; 

• F/W : Filaments/wisps, and; 
Class 2 — Ancillary Classifiers: 
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• R : Ring; 

• E/SO: Elliptical or Lenticular; 

• B : Bar; 

• CL : Cluster, lumpy H il region, knots; 

Four observers (PRH,HK,MJR,KT) inspected the 2x2 kpc postage stamp images in Figure 
[3] and classified each of the Seyfert cores. We did not use the "normal" classifier, because its 
definition could not be independently inferred from MGT98. In practice, we note that galaxies 
that showed regular spiral and dust features in their core morphology were more often classified 
as F/W. Conversely those with more irregular spiral and dust features was classified as DI. These 
classifications are not mutually exclusive, i.e., galaxies could be classified as both DI and F/W. In 
Table [T] all unique visual inspections are provided. 

The majority (91%) of galaxies were identified with dust features. Irregular dust features (DI) 
were observed in 42% (13/31) of Syl and 57% (31/54) of Sy2 AGN. In contrast, 68% (21/31) of Syl, 
and (31/54) 57% of Sy2 host galaxies, showed regular filaments and wispy features (F/W). Thus, 
by visual inspection, we find that Syl host galaxies are more regular in their dust morphologies 
than are Sy2 host galaxies, while Sy2 host galaxies are more chaotic or irregular in their dust 
morphologies than are Syl host galaxies. 

To reduce ambiguity in the classification of regular and irregular dust features in our galaxies, 
and to provide a second confirmation of the MGT98 relationship, we developed an additional system 
specifically for the classification of the core dust morphology of Seyfert galaxies. This classification 
scheme is defined as follows: 

• 1- "Nuclear spiral" — Distribution of features resembles a flocculent or grand-design spiral; 

• 2- "Bar" — A bar-like feature in emission or absorption extends outward from the center of 
the galaxy; 

• 3- "Dust-specific classification" — The previous designations considered all structure. The fol- 
lowing classifications describe only the quality and spatial distribution of what we consider 
to be dust: 

— Group A: 

- s- "Late- type Spiral" — Dust appears distributed in a spiral pattern throughout more 
than 50% of the image. The "inner-arm" regions appear to be clear of any dust; 

- i- "Irregular" — No visually distinguishable pattern can be identified in the spatial 
distribution of dust, i.e., the dust is patchy and irregular in form; 

— Group B: 
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- m-"High Extinction" — Dust features appear to be of high column density. The 
galaxy appears highly extincted. Dust lanes appear to "cut" through the ambient 
stellar light of the galaxy; 

- l-"Low Extinction" — Low contrast dust is present, but is barely discernible from the 
ambient stellar light. 



In Table [21 we provide the classification using this scheme. If possible, galaxies were classified 
using Class 1 and 2, but all galaxies were classified according to their dust structure (Class 3). The 
sub-groups of Class 3 (A&B) were mutually exclusive; e.g., no galaxy could be classified as 'Sis'. 
Galaxies could be classified by a single Group A and one Group B classification simultaneously 
(e.g., '3mi'). If there was a conflict in classifying dust structure amongst the four co-authors, the 
majority classification is listed in Table EJ If no majority was reached after first classification, the 
corresponding author made the final classification without knowledge of the Seyfert class in order 
to prevent any unintentional bias in the measurement of the Malkan relationship. 

The WFPC2 F606W filter we used in this image classification is broad (AA ~4800-7200A) 
and includes the Ha+[NII] line complex. In principle, this line emission could affect our visual 
classification. In practice, the contribution of line flux to the continuum is relatively minor — 
the contribution of the [Nil] doublet to the total flux in this bandpass using the SDSS QSO 



composite spectrum (jvanden Berk et al.ll200ll ) is ^1%, and we estimate the ratio of the equivalent 
widths, EWjj^/EW[jv//], of these lines to be ~ 3:2. Despite the minor contribution to the total 
observed flux in line emission, the photo-ionization of the gas-rich local medium by the central 
engine can produce significa ,nt "hotspots " at th e wavelengths of these atomic lines, which appear 



as structure in the image. ICooke et al.l (|200d ) has studied an example of this photo-ionization 
structure, the spiral-like "S" structure in one of our sample Seyferts hosts (NGC3393; see Figure 
[3|). Though this emission contributes very little to the total flux in the core, the high contrast 
between these bright emitting sources and the local area could lead to "false positive" classifications 
of dust features. Fortunately, few AGN (~7-8 galaxies, see e.g., MARK3, MARK1066, NGC1068, 
NGC3393, NGC4939, & NGC7682 in Figure [3] of the online supplement) show evidence of these 
emitting structures and these highly localized structures are easy to distinguish, in practice, from 
the stellar and dust continuum. 

In conclusion, we confirm that Sy2 host galaxies are significantly more likely to have irregular 
core morphologies: 58% of Sy2 host galaxies were classified as '3i'. In contrast, only 40% of Syl 
host galaxies were classified as '3i'. Furthermore, 39% Sy2 AGN were classified as '3s' in contrast 
to 53% of Syl host galaxies. The results of our visual classification agrees with the observations in 
MGT98. 

Although visual inspection can be an effective means of classifying the morphology of spa- 
tially resolved sub-structure in galaxies, it has its disadvantages as well. Visual inspection is 
time-consuming, it does not provide a quantifiable and independently reproducible measure of the 
irregularity of structures that can be directly compared with the results of similar studies, and 
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it can be highly subjective. Though guidance was provided to the co-authors on how to classify 
varying degrees of dust structure using the Class 3, such classifications are highly subjective and 
conflicts in classification could arise between co-authors. For example, approximately 55% of the 
visual classifications of dust structure (Table [2]) were not unanimous. This discrepancy can be 
largely attributed to the subjective definition of the Class 3 sub-classifications. In each galaxy, the 
co-authors implicitly emphasized certain dust features when making their classification. In many 
galaxies, whether the authors chose to weight the significance of physically small or large-scale dust 
structure could change the structural classification significantly. Consider the case of NGC1365: 
this galaxy was classified with an irregular dust morphology due to the small-scale dust features 
that appear to dominate the visible sub-structure in the core. But, authors who (subconsciously 
or otherwise) emphasized the broad dust "lanes" in the north and (to a lesser extent) south may 
classify the core as having a "spiral" dust morphology. Neither classification is necessarily incorrect 
— the broad dust lanes are clearly associated with the prominent spiral arms in this galaxy when 
viewed in full scale. These complicating factors can weaken any conclusion drawn from the visual 
classification of galaxies. 

In recent decades, as image analysis software and parametric classification techniques have 
become prevalent, the astrophysical community has implemented autorn ated methods for galaxy 



classification (e.g.. lOdewahn et al 



1996 



Conselice 



2003 



Lotz et al.ll2004l ). By relegating the task 



of object classification to automated software and algorithmic batch processing, these methods have 
gained popularity, because they can significantly reduce the time observers m ust spend inspecti ng 
each galaxy, and can provide a reproducible classification for each galaxy (cf. iLisker et al.ll2008l ). 



Therefore, we extend our original test of the Unified Model to include a quantitative assessment 
of the morphological differences between Seyfert galaxies that we identified by visual inspection. We 
present these techniques in f3]and ^ Our use of quantitative parameters to measure morphological 
distinctions between Seyfert core morphologies can provide a more robust test of the MGT98 
relationship by reducing some of the biases implicit in visual inspection. 



4. Traditional Quantitative Morphological Parameters 



A variety of parameters have been defined in the literature to quantify galaxy morphology. 
These parameters are distinguished by their use of a pre-defined functional form — i.e., paramet- 
ric or non-parametric — to e xpress galaxy m orphology. Some popular non-parametric morpholog- 



ical parameters a re "CAS" (IConselice 



and "Gini-M2o" (jAbraham et al 



2003 



20031. "Concentr ation" , "Asymmetry", and "clumpinesS" ) 
Lotz et al.ll2004l . "Gini Coefficient" and M20, the second- 



ord er moment of brigh test 20% of the galaxy pixels). These methods are not without limitations 
(cf . iLisker et al.l |2008| ) , but they each can be useful for assessing galaxy morphology in a user- 
independent and quantitative way. We chose to use these parameters in this research project, 
because the distribution of dust features in the cores of Seyfert galaxies is unlikely to be well- 
described by a single functional form, e.g., the Sersic function that broadly distinguishes between 
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elliptical and spiral galaxy light profiles. 



Conselice et al.l ([20001) provide the following functional definitions of the CAS parameters. 



The concentration index, C, is defined as: 



C = 51n 

r20 



(1) 



where rgo and r2o are the values of the circular radii enclosing 80% and 20% of the total flux. 
The typical r ange in concentration index values measured for galaxies on the Hubble sequence 
is 1 < C < 5 (jConselicel l2004l : iHernandez- Toledo et al.l |2008| ) . Larger values of the concentration 
parameter are measured for galaxies that are more centrally peaked in their light profiles. 

The asymmetry, A, is defined as: 



\Io{a,b) - h{a,b)\ 



A 



a,b=0 



2 E \Io{a,b)\ 

a,b=0 



(2) 



where x and y correspond to the length (in pixels) of the image axes, Iq is the original image 
intensity, and /$ is intensity of pixels in an image that is the original image rotated through an 
angle of $ (we set ^ = 1 80°). Typically, A ranges from (radially symmetric) to 1 (asymmetric). 



see e.g.. 



Consehcel (120031 1. 



dumpiness, S, is defined as: 
S = 10 



E 

a,b=0 



{Io{a,h)-r{a,b))-B{a,b) 
Io{a,b) 



(3) 



where Io{a,b) is the image intensity in pixel (a,b), I'^{a,b) is the pixel intensity in the image 
convolved with a filter of Gaussian widt h a, and B(a,b) is the estimated sky-background for a 
given pixel. Typically, < S < 1 (see e.g.. IConselicell2003l ). and galaxies that appear to be visually 
"clumpier" have higher values of S. 



Abraham et al.l (120031 ) and lLotz et al.l (120041 ) provide the following functional definitions of the 



Gini~M2o parameters. The Gini parameter is defined as : 

n 

^ = n — u 

jn{n — 1) 



n 



(4) 



where / is the mean over all pixel flux values and n is the number of pixels. This parameter 
measures inequality in a population using the ratio of the area between the Lorentz curve, defined 
as: 



(5) 
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and the area under the curve of uniform equahty (= ^ of the total area). Although this param- 
eter was originally developed by economists to study wealth distribution, this parameter can be 
applied to understand the distribution of light in galaxies. If the distribution of light in galaxies 
is sequestered in relatively few bright pixels, the Gini coefficient approximately equals unity. The 
Gini coefficient is approximately equal to zero in galaxies in which the flux associated with each 
pixel is nearly equal amongst all pixels. In other words, the Gini coefficient quantifies how sharply 
peaked, or "delta— function" —like the flux in galaxies is. Note that this parameter can be affected 
by the "sky" surface brightness estimate assumed by the user, which we discuss in Appendix C. 

The M2Q parameter is calculated with respect to the total second-order moment, M^q^, flux 
per pixel, fj, which is defined as: 



M 



tot 



(6) 



such that: 



M 



20 



log 



Mi 



tot 



, while < 0.2ftot, 



(7) 



where Mj is the second-order moment at a pixel j, and (xcU r) are the coord i nates of the central 
pixel . In general, M20 typically ranges between —3 < M^n S, (Lotz et al.ll2004 . 2008 : Holwerda et al 



20 111 ) . If considered jointly with the Gini coefficient. iLotz et al.l (j2004l ) determined that larger values 
of M20 (with correspondingly smaller values of G) are associated with "multiple ULIRG" galaxies, 
and that M20 is a better discriminant of merger signatures in galaxies. 

We measure these five parameters — CAS and Gini-M2o — to quantify distinctions between the 
distribution of light, which underpins the classifications we first made in (©. 



4.1. Case— specific Implementation of Traditional Morphological Parameters 



The authors of CAS and G-M20 (|Conselicell2003l : [Abraham et all 120031 : iLotz et al.ll2004l . re 



spectively) each defined a method to prepare images for analysis that accounts for systematic issues 
(e.g., compensating for bright or saturated cores of the galaxies). This method of image prepara- 
tion and analysis also ensures that the parameters are measured for the galaxy itself, and that the 
contributions from background emission are minimized. In our analysis, we calculate all morpho- 
logical parameters applying a functional form that is consistent with — or identical to — the form 
presented in the literature. However, we caution that our images and specific science goals require 
us to use an algorithm for image preparation and parameter measurement that differs slightly from 
the published methods. In this section, we outline key differences between our data and methods 
we used and those presented in the literature. 

First, we measured these traditional parameters in images of galaxies observed at fundamen- 
tally different spatial resolutions (see §4.2p . All galaxies in our catalog have been observed with 
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HST WFPC2 at a pixel scale of Cf.'W pix"^. In contrast, CAS and Gini-Mao are often measured 
from images obtained with ground- based telescopes that have relatively low spatial resolution in 
comparison with HST. For example. iFrei et al.l (jl996l ) present images obtained with the Lowell 1.1 
and Palomar 1.5 meter telescopes at ~ 2'.'0 resolution at full-width half maximum (FWHM). This 
data set has been used extensively to test the CAS and Gini-Mao parameters' ability to di scriminate 



between the morphological classes and star-formation his tories of nearby galaxies (e.g., IConselice 



2003 



Lotz et al. 



2004 



Hernandez- Toledo et al.ll2006l . l2008l ). The vastly different spatial resolutions 



between images obtained at ground-based observatories and the HST images implies that the pa- 
rameters that we measure are sensitive to features of fundamentally different size— scales. In fact, 
in ground-based images most of the small-scale structure that we identified and used to classify the 
galaxies is undetected. Thus, parameters that are dependent on the pixel-specific flux values 
(e.g., M20), rather than on the average light distribution (e.g., concentration index), may be more 
sensitive to these spatial— resolution differences because at lower resolution fine-scale structure are 
effectively smoothed out. In Appendix A, we quantify the effect of spatial resolution on the specific 
criteria that we use to characterize the structure of dust features in our Seyfert galaxies. 



In IConselicel (j2003l ) and iLotz et al.l (j2004l ). the CAS and Gini-Mao parameters, respectively, 
are all measured in an image that is truncated at the Petrosian radius of the galaxy. The Petrosian 
radius is defined as the radius (vp) at which the ratio of the surface brightness at Vp to the mean 
surface brightness of the galaxy interior to rp equals to a fixed value, r], typically equal to 0.2. A 
Petrosian radius or similar physical constraint is applied to differentiate between galaxy and sky 
pixels, and ensures that the infiuence of background emission is minimized in the calculation of CAS 
and G-Mao- The mean Petrosian radius measured in the r' -filter (Ao=6166A) of Sloan Digital Sky 
Survey (SDSS) Data Release 7 is ~4.3 kpfl Since we are not interested in the contribution of dust 
features located at radii greater than 1 kpc, we do not use images truncated at rp. Furthermore, 
at the mean redshift of our sample, the WFPC2 PC chip field of view is < 2.8 kpc. 

Unlike observations of the entire galaxy, in our core images we can make the reasonable as- 
sumption that most of the flux that we observe in the core arises from sources or features physically 
associated with the galaxy. Not all pixels are sensitive to the flux arising from the galaxy, though, 
and we use the following method to differentiate between the light arising from galaxy sources and 
other extraneous objects or noise. 



• We set all pixels that occur at the edges and the chip gaps between the WFC and PC CCDs 
in the mosaiced images equal to zero. Furthermore, the center of the galaxy is often much 
(10-lOOx) brighter than the rest of the galaxy, likely due to the AGN emission. To avoid such 
extremely bright pixels from biasing the measurement of any of the automated classiflcation 
parameters, we set a high threshold defined as the average of the inner-most 5x5 pixels for 



^Only 28 galaxies in our catalog were observed in SDSS DR7, available online at 'http: //www. sdss . org/dr7' but 
those galaxies common to our survey and SDSS span a range of morphologies and distances, hence we consider the 
measured mean Petrosian radius to be representative for our catalog. 
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each galaxy. We set the pixel values above this threshold equal to zero in the CAS &: G~M2o 
computations. 

If the functional form of a parameter explicitly required a background term, we set this term 
equal to zero. Our analysis is focused on the cores of each galaxy (~1 kpc; or less than 
0.5 XTp), which are significantly brighter, and have high enough surface brightness, that the 
contribution of background objects can be considered to be minimal. We assume that our 
images include only light from the galaxy itself and background emission from the zodiacal 
(foreground) light, which arises from sunlight scattered off of ~100/im dust grains. Fortu- 
nately, from the generally dark HST on-orbit sky, the actual zodiacal sky surface brightness 
is a simple well-known function of ecliptic latitude and longitude We use mea- 

surements of the zodiacal background from WFPC2 archival images presented by Windhorst 
et al. (in prep.), to estimate the emission from this dust in the F606W band. The average 
on-board HST F606W-band zodi acal sky brightness can be found in Table 6.3 of the WFCP2 
Handbook iMcMaster et al.l (120081 ) , but the values in Figure H] give a more accurate mapping 
as a function of i^'^' & b^'^'. The zodiacal background could not be directly calculated from 
the images, because the galaxy core typically over-filled the CCD. We correct for the zodiacal 
foreground emission prior to image analysis in ^4.21 and ^5.1[ For more details, see Appendix 
C. 



To measu re dumpiness, we include d an additional processing step motivated by the algorithm 
de fined inlHambleton et al.l (|201ll ). Prior to calculating the dumpiness parameter as defined 



Conselice 
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^), we first applied a 5x5 pixel boxcar smoothing to the input image. We 



m 

produced the residual map by subtracting the smoothed galaxy image from the original input 
image. The additional smoothed image was generated by applying a boxcar smoothing kernel 
with the one-dimensional size of kernel defined as: 2.0 x ^ x ^, where £ is the dimension of the 
galaxy image in pixels. By design (see f|2]), the linear size of the smoothing kernel is equivalent 
to I or ~0.67 kpc. If we assume that 4 kpc is approximately equal to the Petrosian radius 
for each g alaxy in the samp le, th en this dimension i s coni parable to the smoothing kernel size 
applied in IConselicel (j2003l ) and iHambleton et al.l (j201ll ) . We tested this assumption of an 
average Petrosian radius, and found that using a larger or smaller value (A = ±2kpc) for the 
linear dimension of the kernel has less than ~1% effect on the measurement of dumpiness. 
In this analysis, we also set all pixels within I'.'O of the galaxy center equal to zero. 



In the subsequent analysis, we removed all zero-valued pixels to prevent those pixels from 
affecting the calculation of any of the parameters. 

Though we use identical — or nearly identical — functional definitions of each morphological 
parameter used in the literature, we are analyzing regions of our galaxies at size-scales that are 
significantly different than have been used in previous research. As a result, we cannot assume 
that our parameter measurements are directly comparable to the CAS and G— M20 traditional 



Conselice 


2003: 


Lotz et al. 


2004) 



- 13 - 



parameters that we derived using the above criteria hereafter as C*,A*,S* and G*-M2q, in order 
to distinguish our measurements from the traditional parameters presented in the hterature. 



4.2. Analytical Results and Discussion 

Figure [5] provides three permutations of the measured G*-M2q-C* parameters. In this figure, 
Syl and Sy2 host galaxies are represented in blue and red, respectively. We use this color scheme 
for all figures provided in the online version of this manuscript to distinguish our measurements for 
the two classes of Seyfert galaxies. It is noteworthy that the distribution of each of these parameters 
spans a range that is comparable to the range of the G, M20, and C measured from ground-based 
images at the lower spatial resolution (0.7 < G* < 0.1, -2.5 < M^^ < -0.5, 2.5 < C* < 5.5). 

In Figure [5]^a) we overplot a dashed line that differentiates "normal" galaxies (which reside 
below this line) frorn starb urst galaxies or ULIRGs (i.e.. Ultra Luminous Infrared Galaxies) as 



defined in lLotz et al.l (120041 ). Four of our Seyfert galaxies are measured to be on or above this line: 
NGC1672, NGC4303, NGC4395, NGC7469. The fact that these galaxies reside in this parameter 
space is appropriate, since these four galaxies are considered to be starburst or circum-nuclear 
starburst galaxies in the literature. However, approximately 32% of the Seyfert galaxies were 
identified in the literature as starburst or circum-nuclear starburst galaxies. Hence, we conclude 
that G* and Af|Q do not effectively discriminate between "normal" and starburst galaxies, as 
these parameters are demonstrated to do in the literature. We note that G-M20 are used to 
distinguish starburst and "normal" galaxies when the complete galaxy morphology is considered. 
The morphology of the galaxy on this scale need not necessarily match with the core morphology 
of the galaxies. 

We can consider the relative distribution of the G*-M|q values measured for our AGN. In 
Figure [5{a), we fit a Gaussian function to the G* and M|q distribution and measure the shape, 
centroid, and peak of this function for both Syl and Sy2 AGN to be comparable. The parameters 
of the fitted Gaussian function are provided in Table [3l 

We can draw similar conclusions as above from the distribution of Af|Q — C* and C* — G* 
presented in Figure [5]^b) and (c), respectively. First, it is noteworthy that C* is well-distributed in 
the same parameter space spanned by the conventional concentration index, when it was calculated 
for the entire galaxy at lower spatial resolution. We fit a Gaussian to the C* distribution measured 
for Syl and Sy2 AGN, and measured comparable values for the centroid and FWHM of each 
distribution (see Table [3|) . 

We perform a two-sample Kolmogorov-Smirnov (K-S) for the Syl and Sy2 distributions to 
test whether these distributions are self-similar. The two-sample K-S test can be used to measure 
the likelihood that two empirical distributions were drawn as independent samples from the same 
parent distribution. We use the K-S test here for two reasons, in contrast to more commonly 
measured statistical parameters (e.g., the statistic): 1) the sample size for each distribution is 
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small, which can lead to an incomplete distribution; and 2) a priori we do not know the parent 
distributions from which the empirical distributions were drawn. We use the IDL routine kstwo 
to measure the K~S statistic, d, which equals to the supremum distance between the cumulative 
distribution functions (CDF) of the input distributions, kstwo also reports the probability statistic, 
p, which is the likelihood of measuring the same supremum in a random re-sampling of the parent 
distributions expressed by the empirical distributions. The K-S test cannot provide any insight 
into the parent distribution(s) from which the empirical distributions are drawn, but it can be 
used to test the null hypothesis that the empirical distributions were drawn from the same parent 
distribution. When the K-S statistic is small or the probability is large (p>0.05), the null hypothesis 
cannot be rejected with confidence. 

The results of the K-S test for the Af|Q and G* parameter distributions are provided in Table [31 
These distributions are indistinguishable for both Seyfert classes. However the K-S test measures 
a slightly larger values of (i=0.38 for the distribution of C*, indicating that the CDFs are distinct. 
The associated probability statistic for C* is small (p=0.01). We conclude that the C* distributions 
measured for Syl and Sy2 are significantly different, and thus are likely to be drawn from unique 
independent parent distributions. This could support the morphological distinction between the 
cores of Syl and Sy2 galaxies that was identified by visual inspection in ^ In contrast, if G*- 
M|q are indeed sufficiently robust metrics for distinguishing the distribution of light in the cores 
of these Seyfert galaxies, then the results of the K-S test suggest that these parameters do not 
quantitatively distinguish the galaxy morphologies of Syl and Sy2 AGN. 

We consider the A* (asymmetry) and S* (dumpiness) parameters independently from the G*- 
M|q parameters, because we believe that these parameters are best able to identify — by design — the 
dust features that we found by visual inspection. In Figure [6l^a), we present the distribution of 
A* measured for our AGN. We did not calculate asymmetry for NGC1058, NGC1386, NGC1672, 
NGC3486, NGC4051, NGC4303, NGC4395, and NGC4698, because the WFPC2 images of these 
galaxies included off-chip regions that were set to zero (see ^4.ip . These regions can seriously 
affect our measurements because asymmetry is calculated by differencing a rotated image with the 
original. In Figure [D^b), we present the distribution of A*&:S*. We fit Gaussian functions to the 
distributions of A* and S* measured for Syl and Sy2, respectively. The best-fit Gaussian function 
to each distribution are provided in Table El The Gaussians' parameters measured for Syl and 
Sy2 galaxies appear to be indistinguishable. We confirm this via a two-sample K-S test for the A* 
and S* distributions. The results of this test are presented in Table [3l We conclude from this test 
that both parameters distributions are likely drawn from the same parent distribution. 

The uniformity in the C*,A*,S*, & G*-M2q distributions also suggests that the Ha+[NII] 
emission arising from the photo-ionization of gas (see ^ does not strongly affect the measurement 
of these parameters. Furthermore, if the A* and S* parameters are suitable metrics for quantify- 
ing the morphology of galaxies, then the results of this quantitative analysis do not support the 
correlation between core dust morphology and Seyfert class established by MGT98 and confirmed 
by our visual inspection in ^ 
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In conclusion, four of the five quantitative parameters {A*,S*,G*, and M|q) measured for 
the galaxies do not support the qualitative conclusions developed from visual inspection. The 
distribution of C* may be specific to the class of AGN, which could support the MGT98 relationship, 
but this parameter is the least-suited, by design, to quantify the morphological distinctions that 
supported the morphology-AGN class correlation. We note that this does imply that the MGT98 
relationship is invalid, but that the use of the C*,A*,S*, & G*-M2q parameters to confirm this 
distinction may not be optimal. 



5. Quantitative Morphology with Source Extractor 



The results of the previous analysis may imply that C* , A* , S* , &i G* -M|q parameters are insuf- 
ficient as tools to distinguish the sub-kiloparsec scale features in AGN. In this se ction, we therefore 



develop an additional non-parametric technique that uses Source Extractor (jBertin &: Arnouts 



19961 . hereafter, SExtractor) to measure the distribution of dust features in the cores of AGN host 
galaxies. 

SExtractor is an automated object detection software package that generates photometric ob- 
ject catalogs. This software is widely used for photometry and star/galaxy separation in UV- 
optical-IR images partly due to the software's speed w hen applied to large imag e mosaics. A review 
of the literature returns more than 3000 citations to iBertin &: Arnouts I (119961). with application s 
extending even beyond astrophysics (e.g., medical imaging of tissue cultures bv lTamura et al.ll2010l ). 
The versatility of SExtractor to detect and measure the properties (e.g., aperture photometry) of 
galaxies motivated us to adapt SExtractor for our purposes. In our study, we use SExtractor only 
for object detection, because the algorithm we outline ( §5.ip and apply ( §5.2p may prevent accurate 
photometry. 

SExtractor lias often been used in th e study of nearby, dusty galaxies (see recent work by 



Kacprzak et al. 



2012 



Holwerda et al.ll2012l . for example). This research does not employ SExtrac- 



tor to directly detect and measure the properties of the absorbers. Rather, SExtractor is used to 
derive the photometric properties of galaxies, and these data are coupled with the dust properties 
of the galaxy (e.g., covering fraction). In §5.11 we adapt SExtractor to directly detect dust features 
that are visible to the eye. Thus, our use of SExtractor to outline the characteristics of dust features 
that are fundamentally seen in absorption is a unique application of this software. 



5.1. Technical Implementation to Identify Dust Features 

In this Section, we outline the method we used to train SExtractor to identify "objects" that 
we visually identified as dust features, effectively using this software to mimic the human eye. 

To detect these objects, SExtractor iiisi calculates a local background, and determines whether 
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the flux in each pixel is above a user-defined threshold detect_thresh. All pixels exceeding this 
threshold are grouped with contiguous pixels that also exceed this threshold. When a sufficient 
number (defined by the detectjninarea parameter) of contiguous pixels are found to meet the 
signal threshhold, the pixel group is recorded as an object in the object catalog. Finally, SExtrac- 
tor measures a variety of parameters (e.g., object center, total flux, size, orientation), and constructs 
a segmentation map of detected objects. 

To detect objects corresponding to the visually detected dust features in the cores of our galax- 
ies, it was necessary to flrst train SExtractor using the WFPC2 images of the Seyfert host galaxies. 
Initially, we used the HLA image of each galaxy — appropriately cleaned of defects as detailed in 
^2] — for object detection. After extensive testing, we could not determine a suitable combination 
of the parameters detectjninarea and detect_thresh that would force SExtractor to identify a 
set of comparable objects to the set of dust features that we visually identifled in Sj3l By setting 
detect_thresh low enough that nearly all visually identifled dust features are recovered, too many 
of these features were broken into multiple unique objects. To alleviate the over-segmentation, we 
increased the detectjninarea parameter. In order to recover the majority of the visually identifled 
dust features though, this parameter must be set unfavorably high; dust features were only detected 
when they were included as a component of a much larger, brighter object. 

Direct detection of dust features with SExtractor is difficult. This can be directly attributed 
to the manner in which SExtractor detects objects. SExtractor is designed to detect peaks above 
the local background. In the images, the local background is bright, and not likely to be smooth 
because it arises from the ambient stellar background and not the astronomical/zodiacal sky. Fur- 
thermore, SExtractor can not detect many of the dust features as they are observed in absorption 
with respect to the background. These absorption features may be brighter than the true as- 
trophysical background, but they are still fainter than the local background. Thus, no optimal 
SExtractor parameters can be defined to exclusively select the visually identified dust features. 

We therefore trained SExtractor to identify objects that more closely matched with dust fea- 
tures identified (f|3]) by coupling object detection using SExtractor with the "unsharp-mask" tech- 
nique. The unsharp-mask is a common tool for image analysis, because it enhances features of 
specific spatial scales. In astronomical images, these features correspond to physical objects, such 
as stars, star clusters, and/or dust clouds. To apply this procedure, we flrst convolved the WFPC2 
images with a Gaussian convolution kernel to create a smoothed image. Next, we divided the con- 
volved image by the original image to produce the inverse unsharp-mask image (hereafter, lUMj]]. 
In principle, if we appropriately deflne our convolution kernel such that it enhances these structures 
of speciflc size-scales corresponding to dust features and apply the lUM, those features should now 
be detected as a positive signal above the local background using SExtractor with the appropriate 



"The unsharp-mask image is typicaUy produced by either differencing or dividing the original image by the 
convolved image. When the contrast between the original and the convolved image is small, as it is in our WFPC2 
images, these two different calculations yield similar results. 
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detection parameters. In Figure 6, we provide an illustration of this technique. In Figure 6a&:6b) 
we show the core image of NGC3081 and a surface map of an inter-arm region. We convolved 
the image with a kernel (Figure 6c), and apply the lUM technique to produce Figure 6d. In this 
figure, it is apparent that the dust features in the region of interest have been enhanced by the 
lUM technique. 

To produce the lUM image of each galaxy, we first assumed that giant molecular clouds 
(GMCs) are physically associated with dust features. We generated a Gaussian convolution kernel 
that is specific to each galaxy, the properties of which were motivated by observations of GMCs. To 
produce the appropriate convolution kernel, we first used the galaxy's redshift from NED to define 
a physical pixel scal e (sp\ pixel kpc~^) of the kernel. The linear size scale of GMCs is typically 



less than 100 pc (see lCasoli. Combes. &: Gerinlll984l : iFukui &: Kawamurall2010l ). so we defined the 



FWHM of the kernel equal to i/sp, initially with £=100 pc. We tested a range of size scales, and 
determined that i=80 pc optimally enhanced the sub-kiloparsec scale dust features that we visually 
identified in J}3l We also determined the appropriate linear size of the kernel to be equal to x/W, 
where x is the length on each axis of the image in pixels. 

We determined optimal SExtractor parameters by an iterative process to find the segmentation 
map that most faithfully reproduced the dust features classified in ^ In this process, we fixed the 
5£'a;iractor parameter detectjninarea equal to 90.0/sp for all objects. We required detect_thresh 
for each object pixel to be at least 1.5o" above the local sky-background in the lUM image. Ad- 
ditionally, we determined that the default values for the SExtractor parameters deblendjithresh 
and deblendjnincont equal to 32, and 0.03, respectively, were sufficient for dust feature detection 
in the lUM image. 



We discuss the results of implementing our method using the optimized parameters in ^5.2[ 
The algorithm we have outlined above for the detection of dust features in absorption in images is 
generic. It is not applicable exclusively to our specific scientific interests. Thus, we prepared all 
IDL procedures that we developed to implement this technique for the public. Readers who wish 
to apply this method to other science topics are encouraged to email the corresponding author. 



5.2. Results and Discussion 

In Figure O we presented a four-panel mosaic of 12 galaxies including the WFPC2 galaxy core 
image and its corresponding SExtractor segmentation map. The first of these images was discussed 
in ^ To produce the second image, we reproduced the segmentation images in DS9 using the 
built-in "SLS" color mapE3- This 256-bit "rainbow" color map including black and white allows 
our eyes to better distinguish between different detected objects. However, when the total number 
of detected objects {Nt) is greater than ~40 even this color map is insufficient to distinguish between 



"more details are available online at http://hea-www.harvard.edu/RD/ds9/ 
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all unique neighboring sources. As a result, many unique objects may appear at roughly the same 
color, although these are not necessarily detected as the same object. This limitation of the color 
map does not affect our calculation of Nt. Furthermore, for most galaxies the segmentation maps 
show a number of objects near the edge of the image. Although some of these edge detections may 
be related to real dust features, we excluded these edge detections in the subsequent analysis and 
discussion. 

A comparison of the segmentation map and the galaxy core images suggests that our general 
SExtractor technique is remarkably successful in recovering only those dust features that we identi- 
fied first by visual inspection. Specifically, the dust feature recovery rate using the lUM technique 
is very good for the majority (>95%) of the catalog. For example, bar and spiral arm-like features 
are well-recovered as unique objects (see, e.g., MARK1330). The fidelity of the object detection of 
the spiral arm features is often high enough in these galaxies (see,e .g., NGC3081) that the spiral 
arm features in the image are entirely reproduced in the corresponding segmentation map. 

Galaxies with relatively many dust features — in both regular or chaotic spatial distributions — 
also appear to be faithfully reproduced in their associated segmentation maps. For example, the 
regular structures in NGC1068 and NGC1066 are detected with SExtractor as are the more chaotic 
dust features, as seen in ES0137-634 and ES0323-G77. An interesting result of this lUM analysis 
is that the objects in some galaxies (e.g., NGC1386, NGC1672) are sometimes limited to particular 
quadrants of what appears to be a disk in the original image. The distribution of dust features 
suggests to the eye that this disk (in which the features are embedded) is moderately inclined 
towards the viewer. A discussion of molecular toroid inclination effects are beyond the scope of 
this work, but we will consider this result in future work. We note here that this disk inclination 
was identified first and only by using the SExtractor technique. 

In images where the stellar light profile is exceptionally smooth and few dust features are 
identified by visual inspection, the lUM technique may detect objects that do not strongly corre- 
late with the dust features visually identified in ^ This may represent a limitation of our lUM 
technique. In Figure EJ we included the images of two galaxies (NGC1058 and NGC3608jl!| that 
represent this small fraction (< 5%) of the catalog galaxies. We do not remove these galaxies 
from the subsequent analysis for completeness and to illustrate to the reader instances when our 
SExtractor technique may be limited in its ability to discern visually identified dust features. We 
caution that object detection in these few galaxies using our lUM technique may be more sensitive 
to local pixel-to-pixel noise variations than it is to signal variations arising from dust absorption. 

Furthermore, in those galaxies that appear to have photo-ionization structures (see the 
number and distribution of dust features using the lUM technique still appears to be very well- 
correlated with the visually-identified dust features (see, e.g., NGC3393 in Figure [3|). In some 



"Only four galaxies— MARK348, MARK352, NGC1058, NGC3608— showed any strong distinction between the 
number, size, and spatial distribution of objects detected with SExtractor segmentation map and dust features noted 
by visual inspection in SjS] 
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galaxies, the potential photo-ionization structure appears to be the brightest structure visible in 
the image. It does not appear that this high contrast between photo-ionization structures and the 
visually-identified galactic dust structure has strongly affected the dust feature detection using the 
lUM technique. Variations in the mean signal across these structures could affect the calculation 
of the local sky background with SExtractor, and thus influence dust feature detection in those 
galaxies with possible photo-ionization emission structures. For example, such variations could 
explain the segmentation of what appears as one chaotic dusty region into two approximately equal 
area dust features along the outer edge of the northeastern "spiral-arm" photo-ionization structure 
in NGC3393. 

We also provided in Figure [3] two measurements of the characteristics of the dust structure 
quantified with the lUM technique. We plot the cumulative number of objects for each galaxy 
contained within circular annuli centered at each galaxy's core for a radius Tc, where rc = nxAr and 
Ar=2.0 pixels. Although we present square images of the galaxy, we only calculate the cumulative 
number for annuli with radii less than 1 kpc in the frame to remove edge detections. Using 
the cumulative object number distribution, we calculate a half-object radius (rhaif) defined as 
the radius (in pixels units) of the annulus that contains the inner 50% of the total number of 
detected objects in each galaxy. This value is provided in physical units (parsecs) with measurement 
uncertainties in Table [2l 

In Figure El^a) we plot the distributions of rhai / • We fit a Gaussian function to the distri- 
bution of rhaif for Syl and Sy2 galaxies and provide the parameters of the best-fit functions in 
Table [H There is no apparent distinction in the distribution of half-object radii between Syl and 
Sy2 galaxies. This is confirmed by a two-sample K-S test, the results of which indicate that the 
parent distributions from which the half-object radii distribution were drawn are not likely to be 
unique. 

Figure [3] also included the object surface density distribution (S) measured for the galaxies, 
which we defined as : 



where N is the number of objects contained within annuli of width equal to 10 pixels. We fit a 
linear function to the object surface density function versus radius and measure the best-fit slope 
(a). In Figure [Hb, we provide the distribution of a measured. We fit a Gaussian function to the 
distribution of a as measured for the two classes of AGN, and measure the Gaussian centroids and 
FWHM to be nearly equal (Table |3|) . The similarity in the distributions is confirmed by a two- 
sample K-S test. Thus, the distribution of the object surface density functions does not appear to 
be unique to the class of AGN. 

In Figure [9^, we plot the number of objects (Nt) identified in each galaxy. We measure 
the centroid and half-width half maximum (HWHM) of both the Syl and Sy2 distributions (see 
Tabled]). We fit a Lorentzian function, rather than a Gaussian function, to better account for 
the broad extension from the HWHM peak to large object numbers in the Nt distribution. The 




(8) 
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centroid of the best-fit Lorentzian function of objects equals to ~18 for both classes of AGN, but 
the mean value of the Syl and Sy2 distributions equals to 47 and 35, respectively. Thus, these 
distributions appear to be significantly different. We confirm this result via a two-sample K-S 
test: we measure d=0.35 and p=0.01, and conclude the empirical distributions of Nt measured 
for the Syl and Sy2 populations are likely drawn from independent parent distributions. If the 
objects detected with SExtractor physically correspond to dust features in the galaxies, then we 
may conclude that Sy2 galaxies are, on average, dustier than Syl galaxies. If we remove the four 
galaxies discussed above for which our SExtractor technique did not appear to detect objects that 
are closely associated with the dust features that we identified by visual inspection, though, we 
measure the K-S test probability statistic for the distributions of Nt equal to p=0.06. In this case, 
we can not reject the null hypothesis, and instead are forced to conclude that the distributions of 
Nt measured for Syl and Sy2, respectively, were likely drawn from the same parent population. 

The MGT98 relationship did not consider the number of dust features explicitly, but the 
assignment of relative degrees of dustiness to galaxies implicitly reflects the number of dust features 
that were identified visually. It seems to the authors that is easier to visually classify the dust 
structure as irregular if it contains a relatively large number of dust features, because the eye 
can more readily identify absorption patterns and divergences therefrom. Thus, the mean Nt 
measured for Syl and Sy2 (Figure [9^) may support — albeit indirectly — the MGT98 relationship 
in a quantitative way. 

In FigurelQb and[9h, we also provide the covering fraction (/c) and the average number of pixels 
{Np) associated with objects detected by SExtractor. We fit a Gaussian function to the distributions 
measured for each of these parameters, and observe no distinction between the Gaussian centroid 
or FWHM measured for Syl and Sy2, respectively (see Table H]). We confirm the similarity between 
the measured distributions by a K-S test, and conclude that these distributions are likely drawn 
from the same parent distribution. These results would not support the MGT98 relationship, or at 
least not demand it. 

Throughout this work we have co n sidere d the results of our only in the context of the Unified 



Model, as outlined in lAntonucci et al.l () 19931 ). We restricted our discussion of these results to this 



context, in part, because we were motivated in this work to extend the analysis first presented in 



Malkan. Gorjian and TamI (j 19981 ). in which the authors make a similar assumption on the nature 
of AGN. The assumption of this model is still fair; despite extensive debate (see the Model 
provides a remarkably robust explanation for the observed diversity of AGN^^I. But this model 
is not without rivals. For example, the "clumpy torus" model reduces the thick, dusty torus — 
the inclination of which gives rise to the observed dichotomy of Seyfert-type AGN — to distinct 
individual dust clumps that are generally distributed about the central engine. In this model, the 
AGN type that one observes is not a "binary" function of perspective; rather, the probability of 



'"'"If only because many of the systematic considerations of the Unified Model are still, regrettably, limited by large 
measured uncertainties; cf. Guainazzi et al. 2011 
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observing a Type 1 AGN decreases as the viewer moves towards an "edge-on" perspective but never 
reaches zero. We observe a core region that is hundreds of parsecs beyond the toroid, though. Thus, 
a full interpretation of our results in the context of this model is beyond the scope of this project 
and we reserve that discussion for future work. 

6. Conclusions and Summary 

Recently, multi-wavelength high spatial-resolution images and spectroscopy of AGN-host 
galaxies have revealed that the characteristics (e.g., morphology, dust, properties of the emission 
line regions; see 311) of these galaxies may not be consistent with the Unified Model of AGN. In 
particular, HSTWFFC2 imaging of the cores of local {z < 0.035) Seyfert-host galaxies established 
that the morphologies of dust features in these galaxies correlates with the AGN class (MGT98). 
We investigated this trend by visually inspecting a catalog of archival WFPC2 F606W images of the 
cores of 85 local {z < 0.015) Seyfert galaxies, and classified these AGN by the presence and distri- 
bution of dust features in each. We determined that Sy2 galaxies were more likely to be associated 
with galaxies whose core dust morphology is more irregular and of "later-type" morphology. Thus, 
our visual classification of our Seyfert galaxies confirms the qualitative morphological relationship 
established by MGT98. We concur with the conclusion of MGT98 that — if this morphological 
relationship is indicative of a fundamental distinction between the subclasses of AGN — this result 
weakens the central postulate of the Unified Model of AGN. 

We extended the study of this morphological relationship established through the qualitative 
visual method by re- analyzing the images using quantitative morphological tools. First, in §4.11 
and §4.21 we developed and measured the C*,A*, S* and G*-M2q parameters for the galaxies. The 
distribution of these parameters as measured for Syl and Sy2 AGN did not strongly distinguish 
between the Seyfert class and morphology of the host galaxy. We determined that the parameter 
distributions for Syl and Sy2 AGN are likely drawn from the same parent distribution using a two- 
sample K-S test, with the exception of the concentration C* parameter. In comparison with the 
other four parameters in this set, though, C* is the least sensitive parameter to the morphological 
features identified and classified by visual inspection. We conclude from this analysis that no 
strong morphological distinction exist between the cores of the Syl and Sy2 AGN host galaxies. 
This conclusion conflicts with the established MGT98 relationship which we have visually confirmed 
in g3 

In order to resolve these apparently conflicting conclusions and to address the possibility that 
the C*,A*,S* and G*-M|q are less effective tools for characterizing sub-kiloparsec scale dust 
features, we developed a new method to quantify the core dust morphologies of the AGN galaxies. 
This method combines SExtractor with the lUM technique. We applied this method and found 
that the distributions of the average number of detected dust features in Syl and Sy2 AGN may 
be different. Thus, we have measured a quantitative distinction between Syl and Sy2 AGN that 
supports the MGT98 relationship. Yet, there was no concordance between this result and any other 
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result (i.e., the radial distribution, size and covering fraction of dust features) derived from this 
quantitative method. We therefore cannot strongly distinguish between Syl and Sy2 AGN on the 
basis of their core morphologies using this quantitative method. 

In conclusion, we studied the relation between the host galaxy core (r < 1 kpc) morphology 
and the Seyfert class of AGN. We find evidence to support this conclusion from the results of a 
quantitative assessment of the core dust morphology using existing and new methods. However, 
we can not strongly distinguish between the core dust morphology of Syl and Sy2 AGN using 
the complete set of morphological parameters. Thus, we are reluctant to suggest that the Unified 
Model of AGN must be significantly modified to accommodate the results of this qualitative and 
quantitative analysis. In the future, better and more internally consistent qualitative and quanti- 
tative methods need to be developed to elucidate the true nature of the cores of Seyfert galaxies. 
We expect that JWST will be able to do a similar investigation of this topic at longer wavelengths, 
thereby penetrating deeper into the Seyfert galaxies central dust. 
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Appendix A. Spatial Resolution: Ground vs. Space based imaging 

We implicitly assumed throughout this paper that HST images are necessary to conduct the 
quantitative morphological analyses. If lower spatial resolution ground-based optical images could 
be used instead of the high spatial resolution i^^T images, we could significantly increase the number 
of galaxies we can consider. The SDSS archive, for example, could provide images of hundreds of 
local AGN. 

We downloaded SDSS r' images for 7 AGN that were in both the SDSS DR7 archive and 
our catalog presented in fj2j We made thumbnails of the core (r <1 kpc) SDSS images of each 
galaxies and measured C*,A*, and S* parameters using the same techniques outlined in §4.11 for 
each galaxy. In Figure [TOl we compare these measurements with those presented in which were 
measured in HST F606W images. It is apparent from this comparison that A* and S* cannot 
effectively discriminate between the morphologies of the SDSS galaxies. This result confirms that 
the quantitative morphological analysis we performed above requires the high spatial resolution 
HST images. 
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Appendix B. Size— Scale Relation 

Two galaxies that are identical in every possible way (e.g., morphologically), but at signifi- 
cantly different distances from an observer, will appear different in images obtained with the same 
telescope, because each CCD pixel covers an intrinsically larger physical area in the more distant 
galaxy. As a result, the dust features in the more distant galaxy are less well-resolved spatially. 
Our catalog includes galaxies in the range between 0.001 < z < 0.015, or equivalently a factor of 
10-15 in physical distance. 

We selected six galaxies— NGC1068, NGC3185, NGC3227, NGC3608, NGC4725, NGC4941— 
with morphologies that represented the diversity of morphologies in our catalog. These galaxies are 
all at distances ~ 15 Mpc, and we use these galaxies to quantify the extent to which we are able 
to identify or measure dust features in our catalog galaxies as a function of distance. We do not 
use the nearest galaxies (D < lOMpc) in our catalog because these galaxies include large off-chip 
regions that significantly affect the measurement of asymmetry. 

We then rebinned each of these galaxies to a pixel scale, s, such that : 

„ Dgal N 
s = ^galx— , (1) 

i-'z=0.015 

where I is number of WFPC2 O'.'IO pixels spanning 1000 pc at the physical distance, D, to the 
galaxy and -0^=0.015 corresponds to the distance to a galaxy at the upper redshift range of galaxies 
in our catalog (63 Mpc). 

We measure C*,A*,S* and G*-M|q for these artificially-redshifted galaxies and compare the 
measured values with the original measurements ( ^4.2p . This comparison is presented in Table [5] 
as (5 = ^^y'^^ , where X and Y are the morphological parameters measured in galaxy images at Dgai 
and artificially redshifted to -0^=0.015 • 

In general, the measurement of these parameters does not seem to be strongly affected by the 
relative distance of the galaxy, at least over the relatively small redshift range that we consider 
in this project. For all parameters, 6 is much smaller than the measured dispersion in the range 
of parameters measured in ^4.21 We conclude that range of measured parameters (see Figures 
andE]) are indicative of morphological distinctions between the cores of the sample galaxies, as we 
assumed in the discussion in ^4.2i 



Appendix C. Sensitivity of measurements to the estimated sky background 

Windhorst et al. (in prep.) measured the surface brightness of the zodiacal background as 
function of i^'^- & ft^?- from 6600 WFPC2 F606W and F814W Archival dark-time images. We 
reproduce these measurements for the F606W zodiacal background from Windhorst et al. (in prep.) 
in Figure HI In i j4.1l we estimated the surface brightness of the zodiacal background along the 
line-of-sight to each galaxy in our catalog. We made the reasonable assumption that the only 
background emission present in the images in the cores of our catalog galaxies arises from the 
zodiacal background and then corrected for this background alone in each image. 
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In this section, we measure the uncertainty in our measurements of C*,A*,S* and G*-M2q 
and the object surface density distribution associated with our assumption of the brightness of the 
background in the images. In general, the brightness of the background is determined to have a 
minimal effect on these parameters, with the notable exception of G* , and to a lesser extent, the 
a slope parameters. In Figure 10, we compare the measurements of G* for galaxies corrected for 
a zodiacal background equal to: a) zero (G^); b) the Windhorst et al. background (G^); and c) a 
hypothetical background 10 times larger than the measured in Windhorst et al. The latter estimate 
of the background emission is highly unlikely in any HST image (see Figured]). We assume such a 
large background here only to provide an upper extremum to our measurement of the effect of the 
background surface brightness assumption. 

The dispersion measured for most parameters, i.e., G*,A*,S* and M|q, for different estimates 
of the zodiacal surface brightness was small (<1%). There is a large dispersion between G* and 
G*. We attribute this dispersion to the removal of relatively faint pixels from the measurement of 
G* as increasingly larger values for the sky surface brightnesses are subtracted from the images. 
This has the net effect of (artificially) enhancing the flux associated with relatively higher signal 
pixels, which increases G* . 

We note that there is a modest increase in the measurement of a (the slope of the object 
surface density function) , when comparing cases (a) and (c) . The net effect of background on this 
measurement is less than 5%. Hence, blindly adopting the most likely zodiacal sky-brightness as 
a function of & b^'^' — when this background is not directly measurable — is an acceptable and, 
in our case the only viable, approach. 

Appendix D. lUM Technique: Dust Feature Detection Threshold 

The detection of dust features with Source Extractor is explicitly dependent on the detection 
parameters defined by the user in the configuration file. Here we discuss the typical contrast level 
of the dust features, relative to the "sky background" in the images, which SExtractor detected for 
those parameters outlined in ^5.11 We define the "contrast" as: 

Contrast = ~ ^—^v x 100%, (2) 

I dust ~l~ J meansky 

where fdust is the flux associated with a detected object using the lUM technique and f meansky is 
the average sky value measured in a uniform "sky" region drawn from the core image. 

We measured the contrast parameters for two representative galaxies in the sample, NGC3081 
and NGC3608. The lUM technique appears to work very well in detecting the dust clumps in 
NGC3081, whereas NGC3608 was largely devoid of dust clumps according to the visual inspection. 
For each of these galaxies, we measured the contrast values for three detected dust clumps, using 
two relatively large but smooth "sky" regions (Area~100 — 200 sq. pixels). The mean contrast, 
{fdust = /i the average flux associated with the dust feature) measured for NGC3081 and NGC3608 
equals 6 and 2%, respectively. Assuming fdust equal to the flux of the brightest pixel in each of 
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the dust features, the mean contrast is measured to 12% and 4% for the two galaxies. We measure 
the relative height of the mean flux associated with the dust features above the mean sky equal to 
50-90 XfTsfcy for NGC3608 and NGC3081, respectively. We note here the fainter sources could be 
detected if the SExtractor detection parameters are revised, but this could introduce more "false 
positive" dust feature detections and would fragment coherent visible structure. 
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Figure 1: The distribution of galaxy morphologies compiled from NED. Two Syl AGN were not 
classified in NED. Though these galaxy morphologies are defined for the entire galaxy — not the 
core region which we are investigating — the similarity of these distributions confirms that that 
any distinction that we draw between these classes of AGN is not likely to be attributed to the 
full galaxy morphology. Furthermore, neither class of AGN is biased to a particular class of galaxy, 
nor are we biased generally by our selection criteria towards fundamentally less-dusty galaxy types 
(i.e., early-type galaxies). 
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Figure 2: The FIR flux distribution of the catalog Seyfert galaxies from the IRAS Faint Source 
catalog (available via NED). The Syl distribution has been scaled to match the peak in the Sy2 
distribution and both samples each had one AGN with measured FIR flux greater than 40 Jy 
(illustrated by arrows). We fitted an exponential function (oc exp[— //r]), where r=3.8 &: 6.7 for 
Syl and Sy2, respectively. We did not select Seyfert AGN on the basis of their FIR properties, 
but the samples appear to be generally similar, with the caveat that the sample has a known bias 
towards Sy2 AGN. 
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Figure 3: From left to right, we provide the WFPC2 F606W postage-stamp image of each of the catalog that was used to 
classify galaxy morphology qualititately (Jj3]) and quantitatively (S|4] and Next, we provide the segmentations maps that 
were generated using the inverse unsharp-mask method defined in ^5. H are provided. Finally, we provide the cumulative number 
function of objects measured for radii less than 1 kpc and the half-object radius as well as the object surface density, defined 
as the number of objects per annulus and the best-fit slope a. We discuss each of these data products at length in N5.2I 
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Figure 4: Windhorst et al. (in prep.) measured the surface brightness of the zodiacal background 
as a function of ecUptic latitude and longitude using ~6600 dark-orbit, archival F606W WFPC2 
images. We estimate the surface brightness of the zodiacal background along the line-of-sight to 
our catalog galaxies, and correct for this zodiacal emission by subtracting the background from the 
core image in ^4.11 
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Figure 5: The C*,G*, & Af^Q parameters were defined in ^ Here, we plot the measured pa- 
rameters for Syl and Sy2 AGN as blue circles and red s quares, respective ly. In (a), we overplot 
the empirically-defined line (dashed) that was defined bv lLotz et al.l (j2004l ) to separate "normal" 
galaxies from Ultraluminous Infrared Galaxies (ULIRGs), but find that this line does not strongly 
differentiate starburst-type galaxies from "normal" galaxies. The distributions of each of these 
parameters appear indistinguishable for Syl and Sy2 (see §4.21 for more details). We confirm this 
with a two-sample Kolmogorov-Smirnov test. The measured K-S parameter is large for both the 
G* and M^g distribution (d=0.28 and 0.29, respectively), but the associated probabilities are also 
both large (p=0.09 and 0.08). Therefore, we do not reject the null hypothesis that both distribu- 
tions are drawn from unique parent distribution. However, the K-S test for the distribution of C* 
does suggest that the measured distributions for Syl and Sy2 AGN are drawn from unique parent 
distributions ((i=0.38 and p=0.01). We discuss these results in f|6l 
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Figure 6: The distribution of ^4* and 5*, our non-parametric measure of asymmetry and dumpiness 
(as defined in Q measured for our AGN. We fit each distribution with a Gaussian, and measure 
the centroids and FWHM of these distributions — the distributions appear indistinguishable. The 
results of a K-S test suggests that the two distributions are not independent. 
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Figure 7: A cartoon representation of the inverse unsharp-mask technique (see §5.ip for detecting absorption of stehar hght 
by dust and ckimpy structures along the line-of-sight. In panel (a), the 2kpcx2kpc postage stamp of NGC3081 is provided; in 
this black and white image, black indicates relatively high signal and white indicates low signal. The thick black square in this 
figure emphasizes a spiral arm and inter-arm region with interesting dust features and morphology. A surface map "zoom-in" 
of this region is provided in panel (b); the arm is indicated by the "trough" extending in an arc from east to west in this image. 
To produce the inverse unsharp-mask image, we smoothed image (a) with a representative kernel (panel c) and divided the 
convolved image by the original image. In panel (d), we provide the unsharp-mask surface map of the region in panel (b). It 
is apparent in panel (d) that the spiral arm region where dust absorption was most significant in panel (a) is now sufficiently 
above the background to be detectable using SExtractor defined with an appropriate detection threshold. 
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Figure 8: The distribution of the best-fitting exponential slopes a to the object surface density 
profile and the half object radii of objects or features measured for objects in all Syl and Sy2 galaxies 
detected by SExtractor using the lUM technique (see §5. ID . We fit Gaussian functions to each of the 
distributions, and the results of a K-S test confirms that the parents distribution from which the 
distributions were drawn are likely the same. This suggests that there is no significant difference 
between the azimuthally-averaged spatial distribution of objects, and thus the distribution of dust 
features for the Syl and Sy2 populations appears to be indistinguishable. 
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Figure 9: The relative distributions of three statistics derived from the morphological technique 
discussed in i j5.ll Panel (a) [Top]: The number of morphological features Nt detected in the core 
of each Seyfert galaxy. Panel (b): The distribution of covering fraction fc of dust features in the 
sample as defined in ^5.21 fc is a measure of fraction of the total core image area that is associated 
with detected objects. In general, Syl and Sy2 host-galaxies cover similar fractions of area of 
the host galaxy core. Panel (c): The distribution of the average number of pixels Np (i.e., object 
area= Np x (/.'l^; see ^5.2p . Two galaxies (1 Syl and 1 Sy2) were detected with Np> arrOO, indicated 
by the arrow. In all panels, vertical dotted and solid lines indicate the mean and centroid (measured 
from the best-fitting Gaussian, or Lorentzian in Panel a, function to each distribution) of the 
distributions. The parameters of these fits, as well as the results of the results of the two-sample 
Kolmogorov-Smirnov tests of these distributions, are provided in Table [H Only for the distribution 
of total object number Nt does the K-S test suggest that the empirical distributions were not drawn 
from a common parent population. 
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Q 


.61 


o xj u \^ J 


NGC1566 




04h20m00. 


,4s 


-54d56ml6s 


20.6 


1 


.0 


Q 


.20 


SAB{s)bc 


NGC1667 




04h48m37. 


Is 


-06dl9ml2s 


63.0 


2, 


.0 





.22 


SAB{r)c 


NGC1672 




04h45m42. 


5s 


-59dl4m50s 


18.2 


2, 


.0 


Q 


.16 


SB{r)bc 


IN O Lj ^ 1 1 U 




05h52mll. 


,4s 


0*7™ OOr- 

-U ( dz / mzzs 




2, 


.0 





.23 




M r^* o 1 '3 f; 




07h27m04. 


Is 


1 QA^ T A»n /I T 


30.3 


2, 


.0 





.45 


O A±5 J DC 


i N z u o y 




08h43m38. 


Is 




46 


1 


.9 


Q 


35 




NGC2985 




09h50m22. 


2s 


+72dl6m43s 


18.1 


1 


.9 





.21 


(R')SA{rs)ab 


NGC3081 




09h59m29. 


5s 


-22d49m35s 


32.9 


1 


.9 





.23 


(R)SABO /a{r) 


NGC3185 




10hl7m38. 


,6s 


+ 21d41ml8s 


16.7 


2, 


.0 





.49 


{R)SB{r)a 


NGC3227 




10h23m30. 


,6s 


+ 19d51m54s 


15.8 


1 


.5 





.33 


SAB(s)a pec 


NGC3393 




10h48m23. 


5s 


-25d09m43s 


51.8 


2, 


.0 


Q 


.09 


(R')SB(rs)a? 


NGC3486 




Ilh00m23. 


,9s 


+28d58m30s 


9.34 


2, 


.0 





.26 


SAB{r)c 






Ilh06m47. 


5s 


1 "70^ Q /I A'7r. 


36.5 


1 


.5 





.23 




"\T ri n S3 




Ilhl6m59. 


,0s 


+ ioaUomoOs 


17.2 


2, 


.0 





.18 




M O Q 1 S3 




Ilh32m34. 


,9s 


+ 00uU4mU0s 


13.6 


1 


.0 





.50 


oB{s)a pec 


"\T Q *7 O Q 




Ilh39m01. 


,8s 


-o ( a44niiys 




1 


.0 





.10 


[ix )oti[i Jab 


IN vj o 




Ilh56in28. 


Is 


+55d07m31s 


15.2 


1 


.9 


Q 


.11 


SAB{r)b? 


NGC4051 




12h03m09. 


,6s 


+44d31m53s 


9.61 


1 


.0 


0, 


.25 


SAB(rs)bc 


NGC4117 




12h07m46. 


Is 


+43d07m35s 


12.8 


2, 


.0 


0, 


.63 


so"? 


NGC4303 


MESSIER061 


12h21m54. 


,9s 


+04d28m25s 


21.5 


2, 


.0 


0, 


.10 


SAB{rs)bc 


NGC4378 




12h25ml8. 


Is 


+04d55m31s 


35.2 


2, 


.0 





.06 


{R)SA(s)a 


NGC4395 




12h25m48. 


,9s 


+33d32m49s 


4.37 


1 


.8 


0, 


.16 


SA(s)m? 


NGC4477 




12h30m02. 


2s 


+ 13d38ml2s 


18.6 


2, 


.0 





.07 


SBO(s)? 


NGC4507 




12h35m36. 


,6s 


-39d54m33s 


48.9 


2, 


.0 


0, 


.23 


(R')SAB{rs)b 


NGC4639 




12h42m52. 


,4s 


+ 13dl5m27s 


13.9 


1 


.0 





.32 


SAB{rs)bc 


NGC4698 




12h48m22. 


,9s 


+08d29ml5s 


13.8 


2, 


.0 





.37 


SA(s)ab 


NGC4725 




12h50m26. 


,6s 


+25d30m03s 


16.5 


2, 


.0 





.28 


SAB{r)ab pec 


NGC4939 




13h04ml4. 


,4s 


-10d20m23s 


42.9 


2, 


.0 





.49 


SA(s)bc 


NGC4941 




13h04ml3. 


Is 


-05d33m06s 


15.2 


2 


.0 





.47 


{R)SAB(r)ab? 


AGN-host galaxies (Continued) 



AGN-Host Galaxy, (Continued) 



ID 


Alt. ID 


R.A.(J2000) 


Dccl. {,12000) 


Distance 


Sy. Type 


e 


NED Class 




NGC4968 




13h07m06.0s 


-23d40m37s 


40.8 


2.0 


0.52 


{H,')SABO" 


NGC5135 




13h25m44.lB 


-29d50m01s 


56.8 


2.0 


0.30 


SB(s)ab 


NGC5273 




13h42in08.3s 


+35d39ml5s 


14.6 


1.5 


0.06 


SA0"(s) 


NGC5347 




13h53ml7.8s 


+33d29m27s 


32.1 


2.0 


0.23 


{R')SB{rs)ab 


NGC5427 




14h03m26.1s 


-06d01m51s 


36.1 


2.0 


0.14 


SA{s)c pec 


NGC5643 




14h32m40.7s 


-44dl0m28a 


16.4 


2.0 


0.13 


SAB{rs)c 


NGC5929 




15h26m06.2s 


+41d40ml4s 


34.3 


2.0 


0.21 


Sab? pec 


NGC5953 




15h34m32.4s 


+ 15dllm38s 


27.0 


2.0 


0.14 


SAa? pec 


NGC6221 




16h52m46.1s 


-59dl3m07a 


20.6 


1.0 


0.28 


SB(s)c 


NGC6217 




16h32m39.2s 


+ 78dllmS3s 


18.7 


2.0 


0.16 


(R)SB(rs)bc 


NGC6300 




17hl6m59.5s 


-62d49ml4a 


15.2 


2.0 


0.33 


SB{rs)b 


NGC6814 




19h42m40.6s 


-10dl9m25s 


21.4 


1.5 


0.06 


SAB(rB)bc 


NGC6890 




20hl8ml8.1s 


-44d48m24s 


33.3 


1.9 


0.20 


SA{rs)b 


NGC6951 




20h37ml4.1s 


+66d06m20s 


19.5 


2.0 


0.17 


SAB{rs)bc 


NGC7213 




22h09ml6.3s 


-47dlOmOOs 


24.0 


1.0 


0.09 


SA(s)a? 


NGC7314 




22h35m46.2s 


-26d03m02s 


19.6 


1.9 


0.54 


SAB(rB)bc 


NGC7410 




22h55m00.9s 


-39d39m41s 


24.0 


2.0 


0.69 


SB(s)a 


NGC7469 




23h03ml5.6s 


+08d52m26s 


67.8 


1.0 


0.26 


{R')SAB(rs)a 


NGC7496 




23h09m47.3s 


-43d25m41s 


22.6 


2.0 


0.09 


SB(s)b 


NGC7590 




23hl8m54.8s 


-42dl4m21s 


21.6 


2.0 


0.62 


SA(rs)bc? 


NGC7682 




23h29m03.9s 


+03d32m00s 


71.3 


2.0 


0.08 


SB(r)ab 


NGC7743 




23h44m21.1s 


+09d56m03s 


23.5 


2.0 


0.13 


{R)SBO+(s) 


NGC788 




02h01m06.4s 


-06d48m56s 


56.4 


1.0 


0.26 


SAO/a?(s) 


TOL0109-383 


NGC0424 


01hllm27.6s 


-38d05m00s 


48.7 


2.0 


0.55 


(R)SBO/a?(r) 



Notes- °-: Object ID; : NED preferred object ID; in Mpc using WMAP Year-7 cosmology y<omatsi^^^^^U) . 



Table 2: AGN-Host Galaxy, Catalog 



ID"" 


Morpho. 


Ancillary 


Morpho. 


« (X^ 


) 


Half Object 




Class A 


Class. 


Class B 










Radius 


IR1833-654 


DI,F/W 


B 


2,311 


-0 


45 





21) 


967 


53±54.17 


ES0137-G34 


DI,DC 


CL 


3im 


-1 


51 





23) 


927 


31±37.47 


ES0138-G1 


D-NW 


" 


3im 


— 1 


16 





17) 


612 


37±37.45 


ES0323-G77 


F/W 


R,CL 


1,3s 


-1 


90 





06) 


350 


03±61.09 


ES0362-G18 


D-SW,DC,F/W 


R 


3sm 


-1 


04 





16) 


789 


24±64.02 


ES0362-G8 


DC,DI 


" 


311 


-0 


40 





00) 






ES0373-G29 


DI,D-NE,F/W 


B 


311 


-0 


69 





22) 


823 


03±38.28 


FRL312 


DC, F/W 


B 


2,31m 


-0 


94 





31) 


562 


27±44.66 


FRL51 


DI 


CL 


3il 


-2 


14 





07) 


392 


60±57.75 


IR1249-131 


DC, F/W 


R,CL 


3sl 


-0 


93 





06) 


589 


14±59.56 


IR0450-032 


DC 


R 


3im 


-2 


52 





49) 


604 


99±63.79 


MARK352 


" 


E/SO 


" 


-1 


45 





28) 


296 


90±60.49 


MARK1066 


F/W 


B,CL 


2 , 3sm 


-2 


25 





19) 


688 


47±49.10 


MARK1126 


DI,F/W 


B,R 


2,311 


-1 


14 





32) 


821 


58±43.45 


MARK1157 


DI,F/W,D-NE 


B,R 


2,3sl 


-1 


09 





00) 


533 


15±61.70 


MARK1210 


F/W 


CL,R 


3sl 


-0 


81 





67) 


919 


57±55.01 


MARK1330 


F/W, DC 


R 


3sm 


-1 


53 





24) 


879 


67±36.88 


MARK270 


F/W,D-S,DC 


B 


3sl 


-1 


38 





19) 


474 


55±42.56 


MARKS 


DI,D-NE 


B,CL 


2,31m 


— 1 


36 





26) 


817 


55±55.07 


MARK313 


DI,DC;F/W 


B,CL 


Sim 


-1 


54 





73) 


935 


46±26.98 


MARK348 


F/W 


" 


Ssl 


-1 


68 





38) 


646 


81±61.17 


MARK620 


F/W,DI,D-N, 


B,R,CL 


Ssm 


-1 


46 





35) 


551 


32±25.24 


MARK686 


D-W,DC,F/W 




Ssm 


-2 


40 





42) 


485 


03±57.41 


MARK744 


DI,DC 


R,CL 


Ssl 


-1 


36 





33) 


528 


30±36.61 


MARK766 


DI 


" 


311 


-1 


97 





54) 


727 


96±52.74 


NGC1058 


F/W,DI 


CL 


1,3s 


-0 


90 


1 


28) 


699 


35±7.144 


NGC106S 


F/W,DI 


CL 


1 , Ssm 


-0 


22 





85) 


855 


87±15.64 


NGC1125 


DC,D-SW,DI 


" 


Sim 


-1 


36 





19) 


592 


58±44.69 


NGC1241 


DC, F/W 


CL,R 


Ssm 


-1 


72 





00) 


535 


62±55.09 


NGC135S 


DC,DI 




311 


-1 


62 





24) 


594 


62±54.77 


NGC1365 


DC, F/W 


CL 


Sim 


-1 


23 





94) 


620 


82±22.46 


NGC1386 


D-NW, F/W, DC 


" 


1 , Ssm 


-1 


08 





39) 


543 


61±11.95 


NGC1566 


F/W, DC 


R. 


Ssm 


-1 


00 





22) 


715 


97±20.66 


NGC1667 


F/W, DC 


" 


Ssl 


-2 


34 





46) 


641 


56±61.70 


NGC1672 


F/W, DC 


CL 


l,Ssm 


-1 


68 


2 


35) 


558 


84±18.29 


NGC2110 


F/W,DC,DI,D-N 


" 


l,3sm 


-0 


96 





12) 


566 


33±31.96 


NGC2336 


DI 


E/SO 


3il 


-0 


83 





31) 


794 


94±30.19 


NGC2639 


F/W,DC,D-NE 


B 


Sim 


-1 


34 





02) 


513 


96±45.49 


NGC2985 


F/W, DC 


" 


Ss 


-0 


37 





91) 


718 


04±18.17 


NGC3081 


DI,F/W 


B,R,CL 


1,2, Ssm 


-1 


20 





49) 


738 


19±32.72 


NGC3185 


DC,DI 


R 


31 


-1 


99 





08) 


132 


99±16.73 


NGC3227 


DI,DC 


" 


311 


-0 


87 





27) 


549 


46±15.91 


NGC3393 


F/W,DI,DC 


B,CL 


2, Ssm 


-2 


27 





04) 


472 


74±51.05 


NGC3486 


F/W,DI 


R,CL 


l,Ssl 


-0 


65 


2 


00) 


764 


01±9.387 


NGC3516 


DI 




311 


-1 


55 





04) 


474 


70±36.22 


NGC360S 


" 


E/SO 


" 


-0 


70 





SO) 


851 


63±17.23 


NGC3718 


DI,DC,D-SW 




Sim 


-1 


37 





82) 


492 


53±13.66 


NGC3783 


DI 


E/SO 


311 


— 1 


70 





21) 


610 


91±39.84 


"\T O Q A O O 


T? /\A! "PlT 

r / W .Ui 


R,CL 










69) 


508 


14±15.25 


NGC4051 


DC.DI 




"'"^^^^^ 


-0 


96 





64^ 


800 


02±9.650 


NGC4117 


DI,F/W,DC 


R 


Sim 


-0 


25 


1 


14) 


683 


86±12.85 


NGC4303 


F/W 


R. 


l,Ssm 


-1 


77 


1 


28) 


391 


61±21.50 


NGC437S 


DI 




31 


-1 


70 





17) 


465 


32±35.00 


NGC4395 




CL 




-1 


36 


5 


03) 


526 


02±4.402 


NGC4477 


DC,D-E 




311 


-1 


24 





21) 


540 


95±18.62 


NGC4507 


D-S,DI 




Sim 


-1 


80 





39) 


592 


75±48.20 


NGC4639 


F/W 


B 


2, Ssl 


-0 


71 





58) 


793 


61±14.00 


NGC4698 


DI 


E/SO 


31 


-1 


02 





48) 


664 


53±13.88 


NGC4725 


DI 


E/SO 


311 


-0 


93 


1 


41) 


743 


22±16.58 


NGC4939 


D-W,F/W 


B 


2, Sim 


-1 


15 





16) 


697 


48±42.44 


NGC4941 


D-E,DI 




Sil 


-0 


49 


1 


08) 


863 


44±15.24 


AGN-host gal 


axies (Continued) 



















AGN-Host Galaxy, (Continued) 



ID" 


Morpho. 


Ancillary 


Morpho. 




) 


Half Object 




Class A 


Class. 


Class B 










Radius 




NGC4968 


DC,F/W,D-NE 




Sim 


-1 


35 





09) 


609 


44±40.38 


NGC5135 


DI 


R,CL 


3sm 


-2 


72 





38) 


578 


57±55.80 


NGC5273 


F/W,DC,DI 


_ 


311 


-1 


87 





79) 


212 


56±14.64 


NGC5347 


F/W,DI 


R 


Sim 


-0 


44 





35) 


792 


55±31.96 


NGC5427 


F/W 


R,CL 


l,3sm 


-1 


83 





26) 


789 


46±35.80 


NGC5643 


F/W 


CL 


l,3sl 


-1 


18 


1 


25) 


663 


69±16.48 


NGC5929 


DC.DI 




311 


-1 


06 





09) 


667 


89±34.09 


NGC5953 


F/W 


CL 


l,3sm 


-1 


32 


1 


63) 


667 


67±26.94 


NGC6217 


DI,DC 


CL 


3sm 


-1 


09 





16) 


694 


41±20.58 


NGC6221 


DI,DC,D-SE 


CL 


311 


-0 


78 





21) 


624 


49±18.71 


NGC6300 


DI,D-SW 


CL 


Sim 


-1 


08 





20) 


714 


21±15.25 


NGC6814 


F/W,DC 


B 


2,3sl 


-0 


74 





25) 


585 


77±21.46 


NGC6S90 


F/W 


CL 


l,3sm 


-1 


11 





37) 


767 


87±33.10 


NGC6951 


DI,F/W 


R,B,CL 


2,3sm 


-1 


38 





67) 


545 


79±19.56 


NGC7213 


F/W 




3sl 


-1 


22 





26) 


620 


10±24.02 


NGC7314 


D-E,DC,F/W 




Sim 


-0 


58 





29) 


580 


60±19.62 


NGC7410 


DC,F/W,D-NW 




Sim 


-1 


81 





07) 


450 


93±24.03 


NGC7469 


F/W 


R.B.CL 


3sm 


-2 


75 





07) 


519 


75±66.29 


NGC7496 


DC,DI,D-NW 


CL 


Sim 


-1 


28 





36) 


641 


13±22.63 


NGC7590 


F/W,D-NW 


CL 


Sim 


-1 


28 





27) 


684 


94±21.63 


NGC7682 


DI 


CL 


31 


-2 


25 





32) 


691 


53±69.57 


NGC7743 


F/W 




311 


-0 


68 





34) 


557 


79±23.47 


NGC788 


F/W,D-S 




3sl 


-1 


86 





35) 


415 


99±55.44 


TOL0109-383 


D-SE,F/W 




Sim 


-1 


87 





46) 


668. 


87±4S.05 



Notes- Combined results from qualitative and quantitative analyses are provided. 
Cols.2-3:Morphological classifiers were adopted from ^^aUcai^^Goyiai^^n^^raii^ (^^98) 
and are defined as follows: CLASS A— F/W— Filaments^wTsps^^^^^regina^^ust; 

DC= dust lane passing close to, or bisecting, center; D- [direction] =dust lanes | 
on one side of major axis. ANCILLARY— B=bar; CL=cluster, lumpy H II region, knots; 
E/SO— Elliptical; R— ring. We do not use the "Normal" classifier. For details, [S] 
Col 4.: Class B— Classifiers defined to specifically characterize the , 
dust morphology in the cores of galaxies: 
i: "Spiral"; Absorption"; S:"Bar"; and 

3: "Dust Specific Notes" — s or i: "spiral" or "irregular" &z m or I: "High" 
or "Low Extinction". The full description of the classifiers is provided in Js] 
Col. 5: Slope (q:) of best— fit line to the object surface density profile defined in ^5.2| 

Reduced is provided in parantheses. 
Col. 6:Half-Object Radius (pc) with uncertainty. Here," - ■ ■ " indicates that a galaxy 

had an insufficiently few objects to measure the radius accurately. 
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Table 3: Morphology Comparison I: C*A*S* &: G* -M2Q Technique 



Morpho. Syl Sy2 K-S Test 

Param. Centroid FWHIVI Centroid t'WHM 3 p 



G' 


0.39 


0.26 


0.42 


0.35 


0.28 


0.09 


Mio 


-1.89 


0.74 


-2.01 


1.24 


0.29 


0.08 


3.15 


1.27 


3.57 


2.01 


0.38 


0.01 


A* 


0.08 


0.08 


0.12 


0.09 


0.09 


0.99 


S* 


0.08 


0.13 


0.066 


0.13 


0.19 


0.91 



Notes- The parameters of the best— fit Gaussian function to the distribution 
of morphological parameters measured forSyl and Sy2, respectively, are 
provided hero. For more details and a discussion of the implications of 
these results in our test of the Unified Model, sec 



Table 4: Morphology Comparison 11: S Extractor Technique 



Morpho. 
Param. 



Centroid 



Syl 



Centroid 



Sy2 



a 


-1.27 


0.79 


-1.30 


0.88 


0.14 


0.69 


rhalf 


570 


226 


572 


195 


0.33 


0.43 


Nt 


18 


4 


18 


15 


0.35 


0.01 


U 


0.05 


0.04 


0.04 


0.05 


0.37 


0.16 


iVp 


18 


7 


16 


11 


0.12 


0.98 



Notes- The parameters of the best— fit Gaussian function to the distribution 
of morphological parameters measured forSyl and Sy2, respectively, are 
provided here. The HWHM value is provided for Nt because a Lorentzian, 
Lorentzian, not a Gaussian, function was fit to the measured distribution 
of this parameter. For more details and a discussion of the 
implication of these results in our test of the Unified Model, see ^5.11 
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Figure 10: C*,A*, and S* measured for the core (r<lkpc)of 7 AGN using i?5rF606W and SDSS 
r' images. Line segments connect the measured values for each galaxy, and the vector points away 
from the parameter value measured from the HST image. 

Table 5. Distance Dependence" 



Parameter 


S 


G* 


1.5% 




0.8% 


C* 


10.9% 


A* 


2.2% 


S* 


14.1% 



^ — 5 quantifies the 
dispersion of the pa- 
rameter at native and 
artificially redshift spa- 
tial resolution. See Ap- 
pendix B for details. 
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Figure 11: We discuss in Appendix C the robustness of the parameters to the estimate of sky 
background. G* was measured for the galaxies in images produced for three assumptions of the 
zodiacal background surface brightness equal to : (1) zero, G*; (2) estimated from Windhorst et 
al. (in prep.), Gl, and (3) a (hypothetical) lOx larger than Windhorst et al.,G*. In the left (right) 
panel, we show the measured dispersion (5 = ° ) , where Vindicates measurements in scenarios 
(2) and (3). We measure a significant difference (> 20%) only for scenario (3). See Appendix C for 
more details. 



